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Electronic origin of strain effects on solute stabilities in iron 
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Abstract 
Nonuniform strain fields might induce the segregation of alloying solutes and ultimately lead 
to the mechanical performance degradation of body-centered-cubic (bec) Fe based steels serving 
in extreme environments, which is worthy of investigation. In this paper, two typical volume- 
conserving strains, shear strain (SS) and normal strain (NS), are proposed to investigate the strain 


effects on solute stabilities in bcc iron by first-principles calculations. For solutes in each transition 


metal group, the calculated substitution energy change due to SS exhibits a linear dependence on 
aa the valence d radius of the solutes, and the slope decreases in an exponential manner as a function of 
the absolute difference between the Watson’s electronegativity of iron and the averaged value of each 


transition metal group. This regularity is attributed to the Pauli repulsion between the solutes and 


the nearest neighboring Fe ions modulated by the hybridization of valence d bands, and concluded 
= to be originated from the characteristics of valence d bonding between the transition-metal solutes 
and Fe ions under SS. For main-group and post transition-metal solutes, the considerable drop of 
substitution energy change due to NS is concluded to be originated from the low-energy side shift 
of the widened valence s and p bands of the solutes. Our results indicate that the stabilities of 
substitutional solutes in iron under volume-conserving strain directly correlate with the intrinsic 
properties of the alloying elements, such as the valence d radius and occupancy, having or not 


having valence s and p bands. 
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I. INTRODUCTION 


Body-centered-cubic (bcc) Fe based alloys, e.g., ferritic/martensitic steels, are widely 
used in modern industries, especially in those under special environments like heavy load- 
ing, stress and neutron radiation. As well known, alloying solutes play an important role in 
modifying physical, mechanical and other properties of metals, which certainly show a con- 
centration dependency [1, 2]. For example, experiments indicated that the addition of 2-6% 
Cr to bec iron leads to a dramatic suppression in radiation swelling [3-5], while reduced- 
activation ferritic/martensitic (RAFM) steels with 9% Cr exhibit the least embrittlement 
under irradiation [6]. Great efforts have been made to study radiation induced segregation 
(RIS) of alloying solutes in bcc Fe based alloys, which might considerably change the local 
concentration of alloying solutes. As radiation-produced point defects flow to dislocations, 
grain boundaries, and other sinks, those solutes that preferentially associate with point de- 
fects (e.g., tend to form clusters) or preferentially participate in defect diffusion (due to 
lower barriers) will form coupled fluxes, which leads to buildup or depletion in the vicinity 
of defect sinks [7]. This coupled kinetics, as widely accepted, constitute a basic description 
of the microscopic process of RIS. To date, fairly sound understandings have been reached 
on the RIS mechanisms of some alloying solutes in bcc Fe with the aid of computational 


modeling [8-14]. 


For the convenience of discussion about atomic transport in metallic systems under the 
abovementioned special environments, it is presented here the linear relation between the 
atomic flux J® of a chemical species a and the gradients of chemical potential V uê of all 
species 8 under constant temperature: 


J? = = So Las V pP /(keT). (1) 
B 


Herein, Lag stands for the Onsager matrix, kg for the Boltzmann constant, and T for the 
temperature. As an extension of Fick’s law [15], Eq.(1) was suggested by Garnier et al. for 
describing the flux coupling at work in radiation-induced segregation [16]. For example, in 
the presence of a spatial chemical potential gradient of vacancy, the solute coupled with 
vacancy in kinetics can diffuse under the influence of the nonzero off-diagonal terms in the 
Onsager matrix even if its chemical potential is spacially uniform. Nonetheless, macroscopic 


deformations always exist in bcc Fe based steels serving in stress states, and local strains 


are even more severe in the adjacent regions of extended defects like cracks, surfaces, and 
dislocations. In these cases, nonuniform strain fields may bring about noticeable chemical 
potential gradients of solutes, which motivate solute fluxes under the influence of diagonal 
terms of the Onsager matrix. Hence, for a full view of RIS mechanisms in bcc Fe based 
steels, it is necessary to investigate the effects of strain-induced chemical potential gradients 


of solutes. 


In recent years, there has been some research on atomic kinetics in metals considering 
strain effects. Base on first-principles calculations, Liu et al. demonstrated that the volume- 
change effect on the stability of a substitutional solute in an a-Fe host is scaled by its 
formation volume; and in aluminum when applied with volume-conserving rhombohedral 
strain (RS), the stability of a transition-metal solute directly correlates with two intrinsic 
properties of the solute, i.e., the outermost d radius and the Pauling electronegativity [17, 18]. 
Garnier et al. validated the high accuracy of elastic dipole model in calculating the strain 
effects on migration barriers of solute Si in nickel up to a strain of 1% by comparison with 
first-principles results [19]. Based on classical molecular dynamics simulations, Zhang et 
al. suggested that in bcc Fe matrix the interaction between solute Cu and an intrinsic 
dislocation or loop can be divided into two parts: the elastic interaction via the associated 
stress field and the core effect which presents only within a few lattice constants or a few 
Burgers vectors [20]. Herein, the elastic interaction represents a linear pressure dependence 
of the segregation energy of Cu with the minus size mismatch as the slope, which is actually 
identical to the linear volume dependence of substitution energy in Ref. [17]. However, the 
nature of core effect of dislocations on the segregation energy of Cu might need further 


investigation. 


As far as we know, at present there is a lack of research on solute kinetics in bcc iron 
under the volume-conserving strain. Based on electronic structure calculations, this paper 
will attempt to study on the stabilities of substitutional solutes in bcc iron under two typical 
volume-conserving strains, namely, shear strain (SS) and normal strain (NS). The rest of this 
paper is organized as follows: Section II describes the methodology including the adopted 
strain tensors. Numerical results and discussion are presented in Sec. III. In this section, 
the typicalness of SS and NS is clarified among volume-conserving strains, and then the 


embodied forms of Madelung-like energy change due to SS and due to NS are both deduced 
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for a solute. Based on these groundwork, we make an analysis on the changes of substitution 
energy of solutes respectively due to SS and due to NS. We discuss the mechanisms of how 
SS and NS change the stabilities of solutes at electronic level. An additional discussion 
on the stabilities of some solutes in bcc iron under different strains is also given from the 
viewpoint of nuclear power applications in this section. Finally, a summary is presented in 


Sec. IV. 


FIG. 1: (Color online) Sketches of a bec crystal cell changing into monoclinic (a) and orthorhombic 


(b) shapes when it is applied respectively by shear strain (SS) and by normal strain (NS). 


Il. METHODOLOGY 


All first-principles calculations are performed within spin-polarized density functional 
theory as implemented in the Vienna ab initio simulation package (VASP) [21]. The in- 
teraction between ions and electrons is described by the projector augmented wave (PAW) 
method [22]. Exchange and correlation functions are taken in a form proposed by Perdew 
and Wang (PW91) within the generalized gradient approximation (GGA) [23], while the 
correlation energy interpolation is done by Vosko-Wilk-Nusair method [24]. Binary bcc Fe 
based alloys are modelled by 4 x 4 x 4 bec supercells that each contains 127 Fe atoms and 
one substitutional solute. We use an energy cutoff of 350 eV for the plane-wave expansion of 
wave functions and a I-centered 3 * 3 * 3 k-point mesh for Brillouin zone sampling according 
to Monkhorst-Pack scheme. All atomic relaxation calculations are performed at constant 
volume and shape using the conjugate gradient algorithm. The tetrahedron method with 
Bloch] corrections is employed in all accurate energy calculations. For simplicity, the side 
length of the forementioned bcc supercell in equilibrium is set to be four times of the op- 
timized lattice constant of bec Fe (2.833 A). In the calculation of substitution energy the 


following formula is used: 


Bock = EnFe+1S = x E(n+1)Fe = Es, (2) 


n 
n+l 


where Enre+1s is the total energy of a bec supercell containing n Fe atoms and one substi- 
tutional solute S, E(n+1)re is the total energy of the same bec supercell filled with only Fe 
atoms under the same strain, and Es is the energy per atom of pure crystal S in the stablest 
phase. The two strain modes in Ref.[17], SS and NS, are adopted in this paper. SS and 
NS modes respectively change a cubic cell to be monoclinic and orthorhombic, as shown in 
Fig. 1. The strain tensors are described in Cartesian frame with coordinate axis 1, 2 and 
3 being chosen along [100], [010] and [001] directions of a perfect bec crystal, respectively. 


For the strains in SS mode the following equations exist: eu = €22, €23 = €31 = 0. For 


the strains in NS mode, there exist €12 = €93 = £31 = 0. Tablel lists some components 


of different strain tensors. Noticeably, all these strains are volume-conserving, so that the 
Fermi level and the net charge of ions are unchanged when each of them is applied on a bcc 


Fe12751 (S =Solute) model. 


TABLE I: Components (%) of different strains in SS and NS modes. 


€11, €33, €12 €11, €22, €33 
SS89 0.004, 0, 0.873 NS1 —0.990, 0, 1.0 
SS87 0.034, 0, 2.619 NS3 —2.913, 0, 3.0 
SS85 0.096, 0, 4.370 NS6 —5.660, 0, 6.0 


IHI. RESULTS AND DISCUSSION 
A. Madelung-like energy changes of solutes due to SS and due to NS 


Before further characterization of SS and NS strains, it is necessary and interesting to 
make a force analysis on a solute in the Fej275; system. In this dilute binary alloy system, the 
solute endures two major forces, i.e., the chemical (electronic) interaction and the screened 
Coulomb interaction. Obviously, the solute has chemical interactions mainly with Fe ions in 
the nearest two neighbor shells, but scarcely with those in the third as well as farther neighbor 
shells (>4 A). In metallic systems, screening makes the Coulomb interaction between ions 
short-ranged. According to a study on the relation between Madelung potential Vm and net 


charge qn of ions in Als oLis9 random alloy, an ion has screened Coulomb interactions mainly 


with ions in the nearest two neighbor shells [25]. Herein we present this study in more detail. 
Based on locally self-consistent Green’s-function (LSGF) calculations, Ruban and Skriver 
plotted the Vm versus qn curves for this random alloy by setting the local interaction zone 
(LIZ) respectively including up to the second, third as well as farther neighbor shells, i.e., 
LIZ=3,4 and so on. They reported that the inclusion of more neighbor shells farther than 
the second does not affect the results significantly, and thus LIZ= 3, i.e., including the nearest 
two neighbor shells, may be considered to be converged [25]. Hence, when a strain tensor 
applies on our Fej975, system, the change of En of a solute might be mainly originated 
from the change of interactions with ions in its nearest two neighbor shells. As shown in 
Fig. 1(a), when a bcc Fe crystal cell is applied by SS, one half ions in the first neighbor shell 
are equidistantly closer to the central ion compared with the case in equilibrium state, while 
the other half ions become farther away. When NS applies on a bcc Fe crystal cell, towards 
a corner ion, its second nearest neighbors in [100] and [001] directions respectively become 
closer and farther away, as shown in Fig. 1(b). Table II lists the relative changes of distances 
from the nearest three neighbor shells to a central ion in bec Fe respectively due to SS and 
due to NS. These relative distance changes provide a full view of the change of local lattice 
around a solute influencing the Esup. As shown in Table II, besides the changes of distances 
from ions in the third nearest shell that have slight influence on Esu» of the central solute, 
SS/NS results in changes of the distances from ions in the first /second neighbor shell to the 
solute, but almost no change of the distance from the second/first nearest neighbor shell. 
Hence, SS and NS are typical volume-conserving strains in changing the local lattice around 


a solute. 


TABLE II: Relative change ôd (%) of distances from the nearest three neighbor shells to a central 
atom in bcc Fe due to SS and due to RS strains. The total number of atoms with equivalent ôd in 


each neighbor shell is notified in parentheses. 


SS89 SS85 KEN NS6 
ad: —0.6(4), 0.6(4) —2.8(4), 3.0(4) 0.0(8) 0.2(8) 
Ad: 0.0(6) 0.0(2), 0.2(4) —1.0(2), 0.0(2), 1.0(2) —5.7(2), 0.0(2), 6.0(2) 


dz —0.9(2), 0.0(8), 0.9(2) —4.3(2), 0.1(8), 4.5(2) —0.5(4), 0.0(4), 0.5(4) —2.8(4), 0.3(4), 3.0(4) 


When a bcc Fej75; system is applied by a volume-conserving strain of SS or NS, the 


change of Esu» of a solute obviously satisfies the equation suggested in Ref.[18], 
A Esub = qÈ ARA Sp Abo, (3) 


where o? is the net charge of the solute, A Fiat, is the change of Madelung-like potential, and 
AEc is the energy change corresponding to chemical interactions. For the convenience of 
later analyzing the Abt of solutes due to strain, here we deduce the embodied forms of 
A Flatt respectively due to SS and due to NS. Based on LSGF calculations with atomic-sphere 
approximation, Ruban and Skriver demonstrated that the normalized screening charges (i.e., 
induced charges divided by the net charge of the central ion) at atomic sites in the neighbor 
shells around different substitutional solutes in six different pure metals, follow a single, 
common curve as a function of the distance from the solute in unit of Wigner-Seitz radius of 
the pure metals [25]. Their calculations indicated that in bec Pt the normalized screening 
charges in the first and second nearest neighbor shells around a Cu solute are —0.12 and 
—0.05, respectively. Furthermore, according to their reasoning, in a metallic system the 
screened Coulomb (Madelung-like) potential at the lattice vector R belonging to an ion 


with net charge Q at the origin and the screening charges of Q may be written: 


1 Q(R’) 
Fau = ` l 4 
=E Areo < R- R’ (2) 


where Q(R’) is either the net charge Q at the origin if |R’| = 0 or the screening charge 
at the lattice vector R’ if |R’| A 0 [25]. In our case, a reasonable assumption is that the 
net charges of ions in the same nearest neighbor shell around the solute are equal whether 
the Fe,275; system is under volume-conserving strain or in equilibrium, and we assume that 
they are Qı and Q for the first and second neighbor shells, respectively. Supposing only the 
screening effects of the first nearest neighbor shell are included, the change of Madelung-like 
potential A Fax at solute site between under SS85 and under SS89 is 0.0210Q, — 0.0315Qz2, 
while the Abu at solute site between under NS6 and under NS1 is —0.1072Q; + 0.0358Q2. 
Consequently, we can obtain the embodied forms of Madelung-like energy changes of the 
solute due to SS and due to NS, i.e., (0.0210Q, —0.0315Q2)q° and (—0.1072Q, +0.0358Q2)q° 


eV, where all net charges are in unit of electronic charge e. 
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FIG. 2: (Color online) First-principles substitution energy change A Fup of transition-metal solutes 
in bcc Fe matrix between under SS85 and under SS89 versus the valence d radius ra of the solutes. 
For each transition metal group, A Et, values are represented by the same symbol and the broken 
line is the linear fit to Eq. (6). For clarity, the ideal valence d electron number is notified for each 
transition metal group. Inset: the slope C1 of the linear fit versus the absolute difference Lä dan) 
between the Watson’s electronegativity of Fe and the averaged value of each transition metal group. 


The curve is an exponential fit to data excluding those of Mn- and Cu-groups. 


B. Substitution energy changes of solutes due to SS and due to NS 


We first devote to the discussion on substitution energy changes of transition-metal solutes 
due to SS. Ref.[18] verified a linear relation of substitution energy change AE,,, versus the 
outermost d radius rg and the Pauling electronegativity $2 of transition-metal solutes in 
aluminum under RS, i.e., 


A Esub = Cita + Cop + c3. (5) 


Herein c1, C2 and c3 are parameters. RS brings one half ions of the first nearest neighbor 
shell equidistantly closer to the central solute compared with the case in equilibrium state, 
and the other half ions of this shell farther away in an Al based alloy, and this effect is similar 
to that of SS when it applies on a bcc Fe based alloy. In view of this, we plot the AF sup 
of transition-metal solutes in bcc Fe between under SS85 and under SS89 versus their ra in 


Fig. 2. As can be seen, there is a drastic drop of AE, of Sc- and Ti-group transition-metal 


solutes compared with the case in aluminum between under RS85 and under RS89 [18]. 
Herein the values of rg are the products of the standard Wigner-Seitz radii employed in the 
VASP pseudopotential database and the relative radii of outermost d shells to the Wigner- 
Seitz radii of transition metals in Ref.[26]. This distinct difference, which means obvious 
deviation from the linear relation Eq. (5), might correlate with the tail effect of valence d 
bands of the nearest neighboring Fe ions in Fej27,5, systems, since tail effects of valence s 


and p bands were verified to have little influence on Eq. (5) [18]. 


Herein we present a detailed discussion on the tail charge transfer in inter-transition- 
metal systems, which might be unfamiliar to some readers. Previous studies indicated that 
in metallic systems, the tail charge, originating from the tails of neighboring-site bands, do 
contribute to the net charge of an ion, however, this “charge transfer” is in the opposite 
direction to the traditional Pauling electronegativity notions [27]. Furthermore, it was sug- 
gested that tail charge play an important role in charge transfer between the neighboring 
early transition metal (e.g., Lu, Ti) and late transition metal (e.g., Au, Pt) ions, which is 
obviously owing to d bands [27]. Such d charge transfer was predicted from the experimen- 
tal observations of Au core levels and d bands drop in Au-Ag alloy as well as in some Au 
intermetallic compounds [28]. Watson and Bennett deduced a “first principles” electroneg- 
ativity scale gw for transition metals to describe the propensity to gain or lose d electron, 
which even shows good overall agreement with those well accepted electronegativity scales, 
e.g., the Pauling’s [29]. In our Fej275) systems, the differential charge densities around so- 
lute Lu, Y and Zr in (110) plane clearly exhibit charge accumulation due to the reversed d 
charge transfer, and as an example that around Y is shown in Fig. 3. Herein, the differential 
charge density is defined as the difference between the total charge density of a system and 
the superposition of the isolated atomic charge densities placed at atomic sites [30]. For 
completeness, Fig. 3 also displays three other typical plots of differential charge density in 
(110) plane around solute Hg, Cu and Ta, respectively. According to our observations, the 
differential charge densities around solutes of Sc, Ti, V, Cr, and whole groups of Mn, Fe, Co, 
Ni and Cu are similar to that of Cu, those around Hf, Mo, Nb and W are similar to that of 
Ta, and those around Zn, Cd as well as main-group solutes such like Be and Si are similar 
to that of Hg. Obviously, these differential charge density plots reflect the apparent charge 


transfer between solutes and Fe matrix, which includes the reversed d charge transfer and 
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the intrinsic charge transfer conforming to the Pauling electronegativity notions [18, 29]. 
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FIG. 3: (Color online) Differential charge density on (110) plane for solute Hg (upper left), Cu 


(upper right), Ta (lower left) and Y (lower right) in bec Fe matrix under S889. In all pictures the 


solutes are placed at the center. The unit of charge density is e/A®. 


Based on the above analyses, it is convenient to discuss the influence of d band tails 
from neighboring-sites on the A En of transition-metal solutes in bec iron due to SS. The 
interaction between the on-site d band of a transition-metal ion and the d band tails from 
neighboring-sites was termed “hybridization” by Watson and Bennett [29], which indicates 
its chemical bonding nature. Therefore, in transition-metal alloys, the d band tails from 
neighboring sites simultaneously modify the net charge of the on-site ion and the chemical 
interactions. For our Fej275; (herein, S=transition metals) systems under SS, the two linear 
relations, i.e., gescht and AEg~ra, are not valid, so that Eq. (5) is no longer valid. Actually, 
Watson et al. observed the charge transfer deviating from traditional electronegativity 
notions in dilute transition-metal alloys from Mossbauer isomer shift measurements as early 
as 1981 [31]. Nonetheless, as shown in Fig. 2, for solutes in each transition metal group, the 
first-principles AE.y, between under SS85 and under SS89 can be well fitted to the following 
equation, 


A Esub = Cira + Co, (6) 


where C and C are parameters. Table III lists the fitted values of C1, C2, and the adjusted 
R? denoting the goodness of fits. These results demonstrate that for solutes in each transition 
metal group, the linear rg dependence of A Esu is still valid under SS, which is attributed to 


the Pauli repulsion between the solutes and the closer Fe ions in the first nearest neighbor 


11 


shell, just like the case in Ref.[18]. Nonetheless, a further investigation should be made to 
gain deeper insight into the different slopes of the linear relation between AF, and rq 


corresponding to different transition metal groups. 


TABLE III: Parameters Cı and Ca of the linear fit to Eq. (6) and the corresponding adjusted R? 


values. 


group Sc Ti V Cr Mn Fe Co Ni Cu Zn 


Cı 0.0510 0.0645 0.1110 0.1548 0.1741 0.4805 0.2606 0.1340 0.1307 0.1502 
Cə 0.0092 0.0100 —0.0125 —0.0263 —0.0224 —0.2833 —0.1150 —0.0278 —0.0141 0.0033 
adj.R? 0.951 0.864 0.995 0.998 0.982 0.995 0.870 0.723 0.999 0.942 


Because of equal electronic occupancies in the valence d bands, those transition-metal 
solutes in the same group have similar d band hybridizations with neighboring Fe ions. The 
good fitting by Eq. (6) indicates that the Pauli repulsion energy between the solutes and 
the nearest Fe ions scales linearly with ra for solutes in each transition metal group. Differ- 
ent values of the fitted slope Ci are owing to that the Pauli repulsion is modulated by the 
valence d hybridizations between the transition metal solutes and the Fe ions. If Watson’s 
electronegativity @w correctly describes the valence d hybridization between transition met- 
als [29], Cı might be in a certain correlation with |Aġw], which is the absolute difference 
between @w of Fe and the averaged value of each transition metal group. Apparently, the 
strength of the valence d hybridizations between the transition-metal solutes and the nearest 
Fe ions can be scaled by |Aġw]|. As shown in the inset of Fig. 2, the slopes C4 of the linear 
fits to Eq. (6) except those of Mn- and Cu-groups can be perfectly fitted to the equation: 
Cı = ayexp(—az|Ad@w]|) + a3. The adjusted R? is 0.991, and the fitted values of a1, az and 
a3 are 0.464, 7.323 and 0.051, respectively. The @w scales of Tc and Re (1.92 and 1.95, 
respectively) in dř group are much larger than those of the neighboring Mo and W (1.73 and 
1.8, respectively) in d* groups; meanwhile, those of the d? group (1.88, 1.84 and 1.97 for Cu, 
Ag and Au, respectively) are much smaller than those of the neighboring d? group (2.09, 
2.22 and 2.28 for Ni, Pd and Pt, respectively) [29]. These facts might be the reason for the 
Cı values of Mn- and Cu-groups deviating from the exponential curve in the inset of Fig. 2. 


These understandings might provide an appropriate explanation for the systematic linear 


12 


trends of A Esu» versus rg in Fig.2. Owing to the modulation of valence d hybridizations, 
the fitted lines of from Sc- to Mn-groups are orderly located from bottom up in a stepwise 
manner, and so do those of from Ni- to Zn-groups, as shown in Fig.2. The much steeper 
fitted line of Co-group solutes, resembling that of Fe-group solutes, might be attributed to 
the slight valence d hybridizations of solutes to the neighboring Fe ions. 
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FIG. 4: (Color online) First-principles substitution energy change AF, between under NS6 and 
under NS1 of transition-metal and some main-group solutes in bcc Fe matrix. The data points of 
Cu- and Zn-groups are enclosed in a rectangle of dot-dash line since their AEsub drop in a similar 


manner as those of main-group solutes. 


Now we turn to discuss the AF,y, of solutes in bcc Fe matrix due to NS, as shown in Fig. 4. 
For solutes of Zn-group transition metals and the main-group element Al, Si and As, there is 
a drastic drop of AE between under NS6 and under NS1. Zn-group transition metals, also 
termed as post-transition metals in many literatures, have stable filled outermost d bands, so 
that their charge transfer in bcc Fe matrix confirm to Pauling electronegativity notions just 
like those main-group elements [18]. Magri et al. suggested that in random substitutional 
alloys (RSA), the net charge of an ion is proportional to the number of the nearest neighbors 
of the opposite type, i.e., the so-called correlated charge model (CCM) [32]. Based on CCM, 
Johnson and Pinski developed the “charge-correlated” coherent potential approximation in 
density-functional calculations to incorporate the electrostatic energy for RSA, which gave 
the formation energies very close to the experimental results [33]. Therefore, in our Fej275} 
systems, the net charge of a Fe ion in the first neighboring shell around a solute might 


be rather small, and might be written as, Qı~ — ¢°/8, while Q2~0. The Madelung-like 
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energy change of a solute in the Fej275; system between under NS6 and under NS1 might be 
roughly estimated by ~-+0.0134(q?)?, which scales as (Agp)? (Aġp = $2 — OC) for solutes of 
Zn-group transition metals and main-group elements. However, the AF, of these solutes 
do not exhibit any hint of scaling as (A@p)*, as can seen in Fig. 4. For example, for solute 
Al and Si, the Ab are —0.078 and —0.109 eV, respectively, while the Ad, are —0.22 and 
0.07, respectively. This result indicates that the Madelung-like energy change only plays 
a minor role in AF,,, under NS. Similarly, it might be concluded that the Madelung-like 


2 also plays a minor 


energy change of solutes, approximately estimated by ~+0.0026(q°) 
role in A Esup when SS applies on the systems. Conversely, the energy change corresponding 
to chemical interactions A Eet dominates the AF,» of solutes either under SS or under NS. 
According to Hund’s rule, Cu-group transition metals have stable filled outermost d bands in 
ground states [34], and Mn has a stable half-filled 3d° [35], so that in the first approximation 
they might behave like post-transition metals in bonding in some cases. As shown in Fig. 4, 
for solutes of Cu-group transition metals and Mn, there are also a considerable drop of 
AEsup due to NS, which might be owing to their bonding with Fe ions in a similar way to 
those of post-transition metals. However, the A Et of other transition-metal (except Cu- 


and Zn-groups and Mn) solutes only vary in a much smaller scope when NS applies on the 


systems. 


Here we discuss the electronic origin of substitution energy change for solutes in bcc 
iron when SS or NS applies on the systems. Ruban et al. demonstrated that the crystal 
structure of the host plays a crucial role in the solid solubility via investigations on bec-hep 
solution energy differences of 4d metals in other 4d metals [36]. They stressed that apart 
from the host contribution, the solution energy difference owing to different host structure 
is governed by local effects in the form of the d occupation of the solute bonds. These 
conclusions might guide us to seek the origin of strain effects on solute stabilities in iron, 
since it might be treated as different structures before and after a strain tensor is applied on 
the iron host. Therefore, we calculate the band energies and examine the electronic density 
of states (DOS) for representative solutes (Y, Ti, Cr, Re, Cu, Zn, Si, Bi) and Fe ions of 


the nearest five neighbor shells before and after either SS or NS applies on these alloys. 
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Neighbor Shell Number 
FIG. 5: (Color online) Band energy changes AE pang of Fe ions in the nearest five neighbor shells 
around solute Ti (a), Cr (b) and Cu (c) when the systems are respectively applied by NS6 and 
by SS85. For clarity, only the averaged value per ion is displayed for each neighbor shell, and the 
band energies of all the five neighbor shells are set to be zero when the systems are in equilibrium 


state and denoted by a dash line. 


According to Ref.[37], the band energy of an ion in the alloy is defined as 


deen = >, I En EE, (7) 


where Er and mat E) are Fermi level and density of l state, respectively. Results indicate 
that for all these representative solutes, the band energies of Fe ions in their neighbor shells 
change when either SS or NS applies on the systems, indicating the change of electronic 
structure of the host, which obviously correlates with the substitution energy change of the 
solutes. According to our calculations, for all these solutes, the band energies of Fe ions in 


each neighbor shell drop under SS85 compared with those under zero strain (in equilibrium); 
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the band energies of Fe ions in each neighbor shell drop or increase both in a much smaller 
scope when NS6 applies on the systems for transition-metal solutes except Cu and Zn; 
however, for Cu, Zn and main-group solutes, the band energies of Fe ions in each neighbor 
shell drop almost in the same scope under NS6 as those under SS85. Figure 5 displays the 
band energy changes of the nearest five neighbor shells around solute Ti, Cr and Cu when 


SS85 and NS6 respectively apply on the systems as examples. 


—— d (Ti, SS85) 
——- d (Ti, zero strain) 


s (Si, NS6) 
s (Si, zero strain) 


l-projected DOS (arb. unit) 


— p (Si, NS6) 
-——— p (Si, zero strain) 


-8 -6 -4 -2 0 2 4 
Energy (eV) 
FIG. 6: (Color online) Comparison of [-projected DOS of solute Ti (a) and Si (b, c) when the 
systems are respectively under SS85 and under NS6 with those under zero strain. The scales of 


vertical axes are shown for the convenience of comparing the magnitudes of s, p and d DOS. 


Considering the dominant role of valence d bonding in mixing energy of inter-transition- 
metal alloys [36, 37], we examine the d DOS of transition-metal solutes before and after SS85 


applies on the systems. As the example of Ti shown in Fig. 6(a), the d DOS of all transition- 


16 


metal solutes are considerably dispersed under SS85 compared with those under zero strain. 
This indicates that the valence d bands of transition-metal solutes are widened under S885. 
Just like the case of one-dimensional hydrogen chain discussed by Hoffmann [38], the overlap 
between the valence d bands of the solutes and those of its nearest neighboring Fe ions has 
a positive correlation with the width of valence d bands of the solutes, and the electronic 
occupancy of d bond controls the bonding strength. For solutes in each transition metal 
group, there are equivalent electronic occupancies of the d bonds with Fe ions, which means 
similar bonding strength; meanwhile, as the valence d radius rg of solutes increases, the d 
band overlap under SS becomes larger so that the d band widening becomes stronger, which 
might explain that the AF,,, of solutes increase linearly as a function of rg. For solutes 
in different transition metal groups, the difference in electronic occupancy of d bonds leads 
to different bonding strength. For example, there are 16 electrons in the d bonds between 
Zn-group solutes and Fe ions, leading to the most strong “anti-bonding” (highest energy), 
which accords with that the straight line fitted for the An of Zn-group is located above 
other fitted lines; the d bonds between Sc-group solutes and Fe ions have 7 electrons, leading 
to the most strong bonding (lowest energy), which accords with that the fitted line is located 
at the bottom (see Fig.2). Hence, it might be concluded that the regularity of A Eat of 
transition-metal solutes due to SS originates from the characteristics of valence d bonding 
between the solutes and Fe ions under SS. As the example of Si shown in Fig. 6(b, c), it 
is also found that the valence s and p DOS of solute Cu, Zn and Bi are dispersed under 
NS6 compared with those under zero strain, which indicates the widening of s and p bands. 
In our opinion, the widening of valence s and p bands might be attributed to the changes 
of electronic interactions between the solutes and Fe ions in the second neighboring shell 
when NS6 applies on the systems. Consequently, for Cu- and Zn-group and main-group 
solutes, the widening of valence s and p bands of solutes due to NS6, and exactly speaking, 
the low-energy side shift of the valence s and p bands, results in the considerable drop of 
AE sup of solutes. Therefore, the stabilities of solutes under either SS or NS directly correlate 
with their intrinsic properties, such as valence d radius and occupancy, having or not having 


valence s and p bands. 


At last, on the background of nuclear power applications, we discuss the stabilities of some 


substitutional foreign ions in bcc Fe matrix under strain. RAFM steels were developed from 
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conventional ferritic/martensitic steels via replacing Mo and Nb contents by W and Ta 
[6]. However, under SS and hydrostatic strain, our results indicate that Mo/Nb is a little 
more stable than W/Ta as a substitutional solute in bcc iron, respectively due to smaller 
valence d radius (see Fig. 2) and smaller formation volume [17]. Under NS, solute W/Ta has 
almost the same stability as Mo/Nb in bcc iron, as shown in Fig. 4. Ishii et al. measured 
yield strength, tensile strength, and fatigue strength at 2x10" cycles of twelve bcc Fe based 
binary alloys, and their results indicated that Mo, Nb, Ti and W have very similar solid- 
solution strengthening effects [39]. Noticeably, according to our calculations, solute Ti is 
much more stable than Mo, Nb and W in bcc Fe matrix under SS (see Fig. 2) and under 
hydrostatic strain [17], and has the same stability as Mo, Na and W under NS (see Fig. 4). 
Obviously, two factors lead to the excellent stability of solute Ti in bcc Fe matrix under SS: 
one is the much smaller valence d radius and the other is the much stronger d hybridization 
with Fe ions in the first nearest neighbor shell. Therefore, it might be feasible to prolong 
the service life of bcc Fe based steels serving in stress and radiation conditions by increasing 
Ti content and decreasing other solid-strengthening elements, on the premise of maintaining 


the mechanical strength. 


AE ub(€V) 


Al Si As Pb Bi 
Solutes 


FIG. 7: (Color online) First-principles substitution energy change AF, between under SS85 and 
under SS89 (red triangle) as well as AEsup between under NS6 and under NS1 (black square) of 


some main-group solutes in bcc Fe matrix. 


As structural materials in accelerator driven systems (ADSs), bec Fe based steels suffer 
not only the irradiation of high-energy proton and neutron, but also the dissolution corrosion 
of liquid lead-bismuth eutectic (LBE). In general conditions, Pb and Bi are almost insoluble 


in bcc iron. However, first-principles calculations indicated that they are both aggregated 
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easily around irradiation-induced defects to form X-Pb, and X-Bi,, complexes, where X can 
be helium, vacancy or divacancy [40]. Figure 7 displays the AF, due to SS and due to NS 
for five main-group solutes in bec Fe matrix. The AE, of Pb and Bi between under NS6 
and under NS1 are —0.103 and —0.107 eV, which indicate that just like Al, Si and As, both 
Pb and Bi have greatly increased solubilities in bcc Fe matrix under NS. This result might 
indicate that the dissolution corrosion of LBE is more serious to bcc Fe based steels under 
NS than in equilibrium state. As shown in Fig.7, the AF,,, of Pb and Bi between under 
SS85 and under SS89 are greater than those of Al, Si and As. According to our previous 
study [18], this result might be attributed to that the Pauli repulsion energy, corresponding 
to the overlap integral of 3d bands of the nearest neighboring Fe ions and 5d of the solute 
Pb or Bi under SS, is obviously greater than the overlap integral of 3d bands of Fe and 2p 
of Al/Si, or 3d of As. Additionally, we present the formation volumes of solute As, Pb and 
Bi in bec iron, 1.98, 13.27 and 14.31 A®, respectively, which indicate their different volume 


dependencies of stability under strain. 


IV. SUMMARY 


In this paper, SS and NS are proved to be two typical volume-conserving strains for bec 
metals by force analysis of a solute in a bcc Fe matrix. It is verified that the Madelung-like 
energy change of a solute contributes little to the substitution energy change, while the 
electronic energy change is dominant when either SS or NS applies on the bec Fe matrix. 
When SS applies on the dilute bec Fej27.5, alloy systems, the substitution energy change 
A Esu of solutes in each transition metal group is a linear function of the valence d radius 
rq of the solutes, while the slope decreases exponentially as the absolute difference between 
the Watson’s electronegativity of iron and the averaged value of each group increases. This 
regularity exhibits that the stabilities of transition-metal solutes in bcc iron under SS di- 
rectly correlate with their intrinsic properties, which is attributed to the modulation of the 
valence d hybridization to the Pauli repulsion between the solutes and the nearest neigh- 
bors. It is found that the band energies of the neighboring Fe ions around the solutes change 
when either SS or NS applies on the systems, which correlate with the substitution energy 


change of the solutes. It is concluded that the regular change of substitution energy of 
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transition-metal solutes due to SS originates from the characteristics of valence d bonding 
between the transition-metal solutes and Fe ions under SS. For solutes of Cu- and Zn-group 
transition metals and main-group elements, the considerable drop of Ab due to NS is 


mainly attributed to the low-energy side shift of the widened valence s and p bands. 
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